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ON THE OUTPUT OF ACOUSTICAL SOURCES 


by 

Harold Levine 

Joint Institute for Aeronautics and Acoustics 


§1. Introduction 

The performance of an acoustical source, whose existence and nature 
are revealed by a given inhomogeneous term in a wave equation, depends on its 
environmental aspects; this means, in particular, that changes in the power 
output and directivity pattern for any individual source of prescribed type 
can reflect, in a sensitive manner, those made in the surroundings through 
the disposition of objects, the admission of a background flow or a variation 
in configurational parameters. Isolated source characteristics are relatively 
easy to determine on the basis of linear wave equations and a hypothetical 
uniformity in the setting, although the interaction or coupling of sources 
with their surroundings generally poses a formidable problem. Ingard and Lamb 
(1957) consider the simplest type of interaction problem, with fixed and 
localized sources of either monopolar, dipolar, or quadrupolar nature lying to 
one side of a rigid and indefinitely extended plane reflecting surface; and 
they calculate explicitly the so-called power amplification factors or ratio 
of total radiated power throughout the half-space containing the source to 
that generated by the same source in the whole of a uniform free space. Jacques 
(1971) has subsequently discussed the radiation from fixed mulipolar sources 
situated above an infinite plane surface at which an invariable complex impedance 
condition is chosen; his method of specifying power radiation factors, in common 
with that employed by Ingard and Lamb, rests on the field and energy flux de- 
terminations at great distances from the source. Integrals of the flux over an 
appropriate directional range yield formulae whose evaluation is accomplished by 
numerical processes. 
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Acoustical source problems have come to the fore in contemporary 
developments of aerodynamic sound theory which now encompass both the model- 
ling of the sources themselves and the tracing of effects due to nearby sur- 
faces or structural elements. With their prediction of a substantial rise in 
power output for suitably oriented quadrupolar sources lying close to the 
straight edge of an extended plane, Ffowcs Williams and his collaborators 
(Hall, Crighton, Leppington, and Howe, 1970-) ushered in a more systematic 
study of interaction effects whose objective is, typically, that of estimating 
the wave function or sound level far away from the source. Inasmuch as the 
distant field patterns of virtually all source-surface configurations possess 
a complicated nature it is difficult to obtain analytical measures of the 
total power output through direct angular integration. An alternative to the 
latter procedure is given in what follows, this being arrived at via manipu- 
lation (apart from solution) of the pertinent equations for a well posed linear 
boundary value problem and exemplified by the particular 

Theorem ; Let there be a time-periodic monopole source situated at a 
fixed point P in the homogeneous region exterior to a closed (or infinite) 
rigid source S ; and suppose that the spatial part of the complete wave function 

(p(r) = <j>p(r) + <J>g(r) , (1) 

is realized from a singular component, $p(r) , characteristic of the isolated 
(free) source and a regular component, <f> (r) , attributable to the presence of 
the rigid surface S . Then the ratio, 9* , of the time average power output 
of the source, with and without the reflecting surface, is directly expressible 
in terms of the imaginary part of the regular component at the location of the 
source, viz. 
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^ = 1 + — Im <f> s ( P) . (2) 

The advantage of such a representation for calculational purposes, in comparison 
with a conventional one that involves an integral of the energy flux, or |<f)| , 

over the normal directions to a distant control surface, is readily apparent. 

A first indication of the facility afforded by (2) appears in the problem which 
envisages a point source at distance h from an infinite rigid plane (z = 0) 

+ and the pertinent equations 


P (Source) 


(V 2 +k 2 )*(r)e' i “ t - -o(?-P)e" i “ t , z > 0 


it. 

9z 


0 , z = 0 ; 


S (z=0) * s a s * m P-*- e matter to confirm that 


lk r-P lk r-Q 

, ,**\ e 1 , e 1 x 1 

<J)(r) = H 

47r|r-P| 4ir|r-Q| 


, z > 0 


Q (Image) 


in these circumstances: and utilization 


of (2) furnishes the (familiar) result 


9 > = 1 + -^ Im - 
k 


-> 

ik | P-Q 


4tt | P— Q | 


= 1 + 


sin 2kh 


more expeditiously than does the reduction of an integral 

2 tt/2 „ 

4tt Lim {r / |<J>(r)| *27r sin 0 d0} 

r-*» o 

made precise with the asymptotic form 


4> (r ) 'u 2 cos (kh cos 0) , r -> <» 


of the outgoing wave function. 
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After presenting (in §2) theoretical arguments which justify a 
calculation of the net power output at the source position in particular 
instances, the remaining sections (§3-6) are given over to various more 
complicated applications and extensions having to do with different half- 
plane problems, an impedance boundary condition and, finally, with the effect 
of a background flow in the medium. 

§2. A Theoretical Basis for Local Power Calculation 

Consider a fixed object surrounded by a uniform 
compressible medium extending indefinitely outwards in 
all directions; and suppose that the velocity potential 
Re (<J)(r)e ia)t ) descriptive of acoustical disturbances 
in the medium, associated with a time-periodic and 
isotropic point source at P , can be sought after 
through the inhomogeneous wave equation 

(V 2 +k 2 )<K?) = -6(r-P) . 

Let v(r) and p(r) designate the local velocity and variable pressure, p the 
equilibrium density of the medium, with the expressions 

v = -V(}) , p = -ioopcj) (5) 

in terms of the scalar function ( r ) ; and note the concomitant representation 

for the time-average energy flux vector 

N = ( pv> = Im <f> V<j) , (6) 

where * symbolizes complex conjugation. 


• P 



(4) 
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On multiplying (4) with cf> and rearranging it follows that 


* 2 2 


<T(V"+k")<j) = V.(<fV<fr) - Nl 2 + k 2 \$\ 2 = -6(r-P)<j,*(r) 


whence 


^ Im V • (<p*Vip) = -tflm S(r-P)4>*(r) 


Integration of the latter relation throughout the region V outside the sur- 
face S of the object yields 

P - f Im g / / ^ dS - - f Im / «(?-?)♦*<?) dV 
where 3P is the net power delivered at infinity* In case that either 


or 


then 


according 


3(J>/ 9n g = 0 
<f> = 0 


on S 


JP = - Im / 6 (r-P) {— 


-ik r-P 


V 47T [ r— P ( 

to the resolution (1) , and this implies 


+ <p s (r)} dV 


(7) 


^ Im (f> s (P) 

or 

8P=v/(^) = 1 Im 4» g (P) 

as stated in (2) . 

If the source at P has a dipolar nature, with the unit vector *n 
pointing along its axis, and the inhomogeneous equation (4) is replaced by 

(V 2 +k 2 )<p(r) = -n • V p 6(r-P) (8) 


the counterpart of (7) becomes 
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ff = - ^ Im / n • V p 6(r-?){n • V p - — ~r~~ + <J> *(?)} 


4tt r-P | b 


dV 


= - ^ Im / 6 (r-P) {n * V n • V ~ r~~ + n * V(J> (r)} dV 

V 4tt | r-P | S 


Since 


it follows that 


Lim n • 7 n • V, — k l r ~ F l . i k ! 

?,¥ p |?-?| 3 


and 


TP = 

24tt 


+ ^Im 



g> = 


i + 




(9) 


expresses the dipole power amplification factor. 

For a vertical dipole at P(z = h) in the presence of a uniform rigid 
plane (z = 0) 


-* 3 ik| r-Q| 

<f> q (r) = " aT" , 

S \ 4tt | r-Q I V 


z > 0 


where Q(z = -h) is the image of P relative to the plane; and the application 


of (9) yields 


3P = 


3 d^ . sin k(z+h) v 
k 3 dz 2 2 + h z=h 


= 1 - 



— ^— 7 ) sin 2kh + cos 2kh] 
(2kh) 2 2kh 


in conformity with Ingard and Lamb. The corresponding power representations 
for other multiple sources can be obtained in similar fashion. 

An independent derivation of (2), without recourse to the delta 
function, relies on the detailed versions 
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and 


<t> p (r) 

*S (?> 


ik r-P 
e 

# -> -> 
4tt r-P 


G(r,?) 


- / GOr.r'jK^r') dS' 

S 


( 10 ) 


of the singular and regular parts of the wave function, the latter being formed 
in terms of an as yet unspecified distribution of (secondary) sources on the 
surface of the object, with the density 


Kp(r) = 


3n, 


<K r) 


( 11 ) 


At a considerable distance from S , where 

ikr 


G(r,r*) 6 


-ikn.r 


4irr 


r = r 


-*■ r 

->oo, n = — , r onS 


the total wave function acquires the asymptotic form 

ikr 

4>(r) = <f> p (r) + <P s (r) - A(n) , r ■+ oo 

with an angular or directivity factor 

-> ->■ ->. 

Aft\ - -ikn*P . -ikn»r ,-k 
A( n) = e - Je Kp(r) dS ; 

S r 

and it follows, having regard for the energy flux vector ( 4 ), that 


( 12 ) 


( 13 ) 


* m h s > A(n) l 2 

i n M •* 

~ -xkn»P - -ikn»r „ . 

- 4 ^: / il-2 Re e /e Kp(r) dS + 

S 

+ | / e -ikn Kp(r) dS | 2 } d 

s n 

where df}^ designates an element of solid angle about the direction of n and 
n 

/ dfi^ = 47T . 
n 
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To simplify (14) , and thereby confirm the central place of the 
power formula (2), information about the function Kp(r) is relevant; 
thus, if a (soft) boundary condition, <f> - 0 , holds on S the function 
Kp(r) satisfies an integral equation of the first kind, 


<f> p (r) 


/ G(r,r')K p (r’) dS' 
S 


i c ^-^r ? ' i v?’) + 1 

S 4ir|r-r'| F 


f sin k | r-r 1 
S 4Tr|r-r' | 


r on S 


K^? 1 ) dS’ . 


(15) 


^ * 

The consequence of multiplying in (15) by the complex conjugate function Kp (r) , 
integrating over S and extracting the imaginary part therefrom, is that 


Im S &/<?) ^(r) dS - / £ is. ilH l l ds ds . 

S S 4ir|r-r' | 


~~2 f I 1 e" lkn ’ r K (r) dS | 2 dfi (16) 

(4t r) 2 S ^ n 


inasmuch as 


sin kr 1 r ikn«r ^ 

— k7~ ■ 4? 1 e dn * • 

n 


(17) 


Use of the latter result along with prior expressions for (f> , <}> and the 

i j 

soft boundary condition on S implies that 

* -> -> -v -ik|r-P| 

Im /K (r)(p (r) dS = Im(/K_(r) ^ — - dS) 

S T P S ^ 4tt I r-P I 


! m(/K_(?) dS 

S ^ 4tt| r-P| 


- i /, (?) sis. ilMl ds) * 

S ^ 4tt | r-P | 


= Im(-d, s (P) - 2i /Kp(r) dS)* 


4ir I r-P I 


= Im <f> (P) + Im 2i — — / e “ ikn * P f K p *(r)e ikn ‘ r dS dU 

(4tt)^ T n 

= Im <J> S (P) + Re / d^ e ikn * P / e _ikn ’ r K (?) dS 

(4t r) n S 


( 18 ) 
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and, after (16), (18) are combined with (14), the relation (2) is directly 
forthcoming. 

When line rather than point sources are contemplated, the two- 
dimensional analogue of (2), 

&• 1 + A Im <() s (?) (19) 

is applicable. 


§3. Source Radiation in the Presence of a Half-plane 

As the first of several illustrative problems selected to underline 
the relative ease of power calculations made at the source itself, consider 
the output from a periodic and fixed monopole which interacts with a rigid 
half~plane. Suppose that the straight edge of the half-plane extends along 
the z-axis of a cylindrical coordinate system 
and that 


6 = 0 , 2tt 

on the respective faces, where 0 < r < 00 . 
If the source is located at the point 
P( r O»®o’ 0 ^ the wave or Green's function 
specified by the differential equation 

and the boundary condition 


* 0 / 


« r -V 


Half 


Plane 


6(e-e 0 )<s(z) 




0 = 0 , 2tt 


admits the representation 
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4>(r, e, z ;P) = G(r,0,z;r o , e Q ,0) 


i f ikh " 

1 h + *£ / 

2 I 4ttR_ 4tt q j 


2/rr 0 cos — — /H 1 <1) <k/R_ 2 +5^)' 


ikh + 

+ £ + i^ y 

4irR + 4 t it 


e+e n 

^ cos —/h, <» 


iU 


V+s 2 / J 


where 


R__ 2 = r 2 + r^ 2 - 2rr^ c °s(S+0q) + z 2 


and ^(?) designates a first order Hankel function of the first kind. 


The power factor 8P of this source is, in accordance with (2), 


found from 


1 + X Lim Im{G(r,e,z;r 0 ,e 0 ,0) 


and the outcome 

*1.1 sin(2kr o sl ” V , 1 . 2kr ° J i<« , r , i r 2kr ° COS%J /' 2kr o sin V 2+ « 2) 

- 2 2 2kr sin 6 + 2 / 5 d5 + 2 ' r ■ , ~ J{ 

/(2kr Q sin 0 Q ) Z +? Z 


involves the first order Bessel function J ]L (C) • It follows from (21) along 


with the result 


°° J,(ai4 2 +5 2 ) 

0 iV 


— i \cxru tl / 

r 1 sin ab ^ „ 

/ — at, = ; , a > 0 

, n , — o ^ ab 
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9 > 


^ sln(2krQ sin 0g) 2 00 J 1 (^2kr Q sin 0 Q ) 2 +£ 2 ) 

- 1 + "7 oVv o *■ V f ' ' d£ + 


2 2kr Q sin 0 Q 


A 


2kr Q sin e 0 ) 2 +£ 2 


d? + 0((kr 0 )" 3/2 ) 


= 1 + 


sin 2kh 
2kh 


as kx 


0 


where 


kr^ cos 0g ■+ oo , 0 < 0g < tt/2 


h = r 0 sin 0 Q 


specifies the finite normal distance between the source and the plane; agree- 
ment of the latter estimate for a source far removed from the edge of the half- 
plane with the previous determination (3) for the same source in the presence of 
a full rigid plane is only to be expected. 

An infinitesimally thin rigid screen does not, evidently, exercise any 
influence on the output of a coplanar monopole point source; and the requisite 
value 1 is assumed by (21) when 0^ = tt , inasmuch as 


1 " 2kr o^(0 

2 o 5 


! 2kr 0 J AO 
- - J f -V" d£ 
0 * 


and 


sin (2kr Q sin 0 Q ) 
2kr Q sin 0 Q ” 


->1 , 0 Q ->■ IT . 


The power factor that replaces (21) in the case of a soft half-plane, 
at which c|> itself vanishes, is 
sin(2kr sin 0_) 


<?>=!_ A 


■o 


"O' 


2 2 2k r Q sin 0^ 

1 2 , kr 0 J l ( ® 


2kr rt cos 0 


+ 4 f S 


0 


0 


/(2kr sin 0 n ) 2 +? 2 


d? (22) 


(featuring opposite signs of the second and fourth terms relative to those in 
(21)) and the particular deductions therefrom, namely 
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- 1 - 


sin 2kh 


2kh 


, kr Q cos 0 Q -> oo , h = r n sin 0 n , 0 < 6 n < ir/2 


0 


and 


2kr 


0 J.(5) 

/ — | — d£ , 0Q = tt , 0 < r Q < 


( 23 ) 


are noteworthy. The latter function of a dimensionless argument a = kr^ , 
where r^ is the least distance from the source to the edge of a coplanar 
screen, possesses a regular nature for all a and has the limiting behaviors 
1 , a ■+ a> , a , a -* 0 . 

In two dimensions the source or Green's function which satisfy (4) take 
the respective forms 


<t>(r,6;P) = 


G(r,0;r o ,0 o ) 


= 1 Ji H 

2 14 H 0 


(1) (kR ) + ~- / 

u - 2 tt o 

0+0 

2/cxZ cos — r 

±4 H 0 (k V±£ ' 

0 


0- e O 

2v ^0 COS ~~T ~ /exp [ik/R 


iA 2 +£ 2 ] N 


AW 




f exp [iky^" 


2 +? 2 ]' 


4 ?? 




(24) 


when rigid or soft boundary conditions hold at the trace 0 = 0(2ir) , 0 < r < °° 
of a half-plane; and the designations 


R_ 2 = r 2 + r Q 2 - 2rr Q cos (0 + 0 Q ) 

refer to the (squared) distances between the source, and its image in the plane, 
from an arbitrarily located point r,0 . 

The power factor which bears direct correspondence with (23) , that is 
for an isotropic line source located in the extension, and parallel to the edge, 
of a half-plane proves to be 
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9 = - / 
" 0 


^ ? Si(2kr 0 ) > 


where Si(^) denotes the sine integral, and this also evidences a regular 
behavior, as a function of kr Q , on the range 0 £ r Q < 00 . 

A rather different behavior of the output generally obtains if higher 
order primary sources (with more intense near fields) are chosen; and, spe- 
cifically, this means an indefinite rise in power amplification as the source 
moves closer to the edge of the plane. Consider, for example, a dipole source 
which is oriented normally to the line that issues from the trace of a rigid 
half-plane and passes through its center. 

The appropriate (two-dimensional) power ^ s' 

formula in this circumstance, ^ / 

xfe 6=0 

<a> , . 8 3 , N r 0 2lT 

9 1 + 72 — Im -30 ^S (r 0’ 7T) » 

k r o 

where ^ denotes the regular (secondary) part of a source function defined 
by the equations 


(V +k )G = - 


«<r-r n ) 


_0_ _3_ 

2 30 


6(0+ir) 


= 0 , 0 = 0,2tt , r > 0 


yields, on employing the (monopole) source representation (24), 


1 + i if? (- I H 0 C1 ) OcR_) + f H 0 (1) (kR + ) 

'• kr 0 ; U ^ 0-0Q 

2r 0 cos — 7 exp(ik ^ 2 +C 2 \ 

+ -T-S == 


AW 


0+0Q 

1 2r 0 cos 2 ~~ ( exp (ikv4 + 2 +£ 2 ) 




v4. 2 +£ 2 


0=0 o =TT 
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= 1 + 


8 3 2 fl _ ... 

2 „ n 2 14 J 0 (2kr 0 COS 2 


(kr Q )^ 30 . 


2r sin -r ,, /I 0\ 2 ." 2 ■>. 

g 1 0 2 sxn(k/(2r 0 cos j) +£ ^ 

) — — - / ~ ~ ~ -1 ' ? 

0,2 2 Je-ir 


/(2r Q cos -j)“+r 


= 1 + 


(k r Q ) 


2 30 


3 f kr 0 . 0 T ... .0, 1 0 . 

sm j J 1 (2kr Q sm j) - cos sm 2kr Q 


V 


0 


, 1 2 . Q , 

+ 2? r 0 Sin 9 Q f 


2r 0 sin 2 


cos(k/(2rQ cos y) 2 +^ 2 

7^ 0s2 . _2 

(2r Q cos j) +? 


sin(k/(2rQ cos 

I 0,2 2.3/2 

[ (2r Q cos -j) +5 ] 


sin 2kr 


2kr, 


0 , 4 


"(kr 0 )' 


, n r ,sin T cos T\ j 

+ - / (—3 -J-) dT . 


d? 


J 0=TT 


Since 


2kr o . 2kr o 

/ — k " 0T = ! sin Td( 

e t e 2 t 

. oi 2kr 

sm 2kr rt . .. 0 

0 sm e l - cos t 

2 + ? + 2 ^ t 

2(2kr Q ) 2e e 


dT 


it turns out, finally, that 
sin 2kr 


2kr r 


^ = 


0 , 2 _ . sin £ r 
— j + — Lim J 

27r(krg) 77 £-K) I e e 


sasjiat) 

X 1 

/ 


„ sin 2kr n cos 2kr n 

« ^ Si(2kr_) + £ + — - 

TT O' n _„ x 2 TTkr, 


27r(kr 0 )‘ 


'0 


(26) 


Thus the output factor approaches unity when the dipole source recedes far away 
from the edge (kr^ °°) and, as the estimate 


& = — r- — , kr -> 0 

irkr^ 0 


makes plain, becomes arbitrarily large when the source nears the edge. 
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Before concluding this section an idea may be given of the intricate 
analysis which is in prospect for total power calculations along conventional 
lines. As regards, in particular, the relatively simple problem of an iso- 
tropic line source symmetrically placed with respect to a half-plane, whose 
power factor involves the sine integral [cf. (25)], let it be observed that 
the multiple integral characterization 

„ 2tt '"“O 2.2 12 

0 


f 

IT 0 


/2kr Q sin j 2 
/ e^ dy 


de 


(27) 


is consequent to joint use of the far field estimate 


* 


. ik(r+r cos 6) Sin 2 .2 

i 0 * iu 

-pZT e S e M dy , 

✓5kr 0 


kr » 1 


and the energy flux vector (4) . Taking account of the Fourier series development 

ikr Q cosG sin 2 2 co -Tri(n+|)/2 

e f e dy = e Ee (-1) sin(n+i)9J , (kr n ) 

0 n=0 2 n+| 0 

the power integral (27) goes over to a single sum formula 

00 


*’- 2 E „ J .i< k V 


n“0 n+j 

which, after transformations based on Bessel function properties and a change 

of integration order, reverts to a single definite integral, namely 

u/2 tt/2 f - 2kr 0 COS ^ 

? / J 2n+l (2kr 0 cos d * = 7 / 1 I J J 0 (C) ^ 


2kr 0 cos _1 (^/2kr 0 ) 2kr Q 

«V5) / d *-£/ J 0 <5) cos ' 1 


d? . 
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Since the function 


a 


1(a) = / J (O cos" 1 £ d£ 
0 a 


has the first derivative 
a 


tt/2 


dl £ 

jT = / J n (D ~ " " " / sin \p J n (a sin ij>) dij/ 

0 0 0 


sin a 
a 


and the initial value 1(0) = 0 , the determination 

1(a) * f d 5 = Si (a) 

0 * 

follows, and this secures the anticipated result 


9 - i S1<2kr 0 ) • 


§4, Radiation from a Line Source Near an Edge at which a Kutta Condition Holds 
It is natural to speculate on the link between enhanced sound radiation 
for (multipolar) sources near a sharp edge and the presumption of locally in- 
compressible flow; Jones (1972) undertook, in this regard, to analyze the effect 
on sound excitation by a time-periodic line source near the edge of a rigid 
half-plane when a Kutta condition suppresses the unbounded nature of the edge 
velocity that is characteristic of incompressible and inviscid flow. His model 
problem, wherein the trace of a rigid plane y 

occupies the half-line y ■ 0 , x < 0 , 
envisages both a course at P(x Q ,y 0 ) and the 
existence of a wake (or discontinuity line) 
elsewhere along the x-axis. 


* p(x o ,y o ) 




WAKE 
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The postulated equations for the coordinate factor <J>(r,P) of a 
velocity potential comprise 

(V 2 +k 2 )<(>(r,P) = -6(x-x 0 )6(y-y 0 ) 

, y = 0 , x < 0 (28) 

and 

<f>(x,0+;P) - <f)(x,0-;P) = Ae lKX , x > 0 


wherein the latter (or wake) condition involves two constant parameters A 
and k ; only one of these constants, namely A , can be fixed by means of the 
foregoing equations along with the further hypothesis that 3cf)/ 3y (i.e., the 
normal velocity) is continuous at the wake and the Kutta requirement. The other 
parameter, k , is assigned a positive (though indeterminate) value, as befits 
an outgoing wave or radiation condition* 

Let the superscripts + differentiate the wave factors in y > 0 , 
respectively, and define the function 


vw - |i 


y=0 


> x > 0 ; 


(29) 


then, if y^ > 0 , say, the representations 

4> (+) (x,y;P) = | Hj 1) (k^x-x 0 ) 2 +(y-y 0 ) 2 ) + j H^ 1} (k/(x-x Q ) 2 +(yfy 0 ) 2 ) - 

- j / (k*/(x-x , ) 2 +y 2 )V(x') dx' 


and 


<f> (-) (x,y;P) = | / H^^k/cx-x'jVjVCx') dx' 


(30) 


comply with the first two equations of the set (28) and jointly provide, on im- 
posing the wake condition, a linear integral equation for V , viz. 
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00 ^ 

J H 0 ^ k l x "" xl |) v ( x? ) <ix f = Y Hq 1 ^ (ky/U“X Q ) 2 +y 0 2 ) + lAe ±KX , x > 0 . (31) 

A formal resolution of the boundary value problem is thus achieved by solving 
the integral equation and employing the representations (30). There is no need, 
in fact, to investigate that part of the solution which obtains in the absence 
of a wake discontinuity (A = 0) , since the corresponding source or Green's 
function of a rigid half-plane is known. Application of the Wiener-Hopf method 
to the integral equation that contains the wake term alone, 


f H^(k|x-x' l)V(x') dx' = iAe llOC , x > 0 


iKX 


(32) 


yields, in sequential order, the complex Fourier transform of V(x) , name ly 


v<?> - / e i5x V(x) dx = - f ,£+k - — , c > 0 


iE. 


and thence the function itself, 


(33) 


V(x) = ~ / e" 1?x V(0 d£ - - v£+k f y ^ jSEL 

/7r _oo 4-rr ^ £ + k + ie 


dC 


(34) 


On referring to an alternative form of the latter. 


v( X ) = - ~ ; e 


-iCx 

A+?; 


00 -i£x , 

d? - (K-k) / — 


- 00 A+c (c+K+ie) 


and the result 

— i£x , 00 _ixx 


; £^i . 2e -W2 / dT . 2 /g „ 

v4c+C k /r^k x 

A 

the singular part of V(x) at x - 0 is isolated, viz. 


V(x) = - (K+k) e i(kx-7T/4) , x -> 0+ . 


(35) 



-19- 


A calculation made with the Green’s function of a rigid half-plane 
lying on the negative x-axis reveals that 


v(x)= 2 kiel ' 


2^ 0 exp(ik»4 2 +r 0 2 -2xr 0 cos (0-0 Q )+5 2 


*^ 2 +r 0 2 -2xr Q cos ce-e > + ? 


dZ 


0=0 


r— a ik(x+r) 
I A 0 O e 0 

2tt / x Sin 2 x + r„ 


6, 


- ^ lkr Q sin 8 0 


2>/Xr 0 C ° S 2 exp(ikv4 2 +r 2 -2xr n cos 6 + £ 2 ) 

! — j — r — 2 2 — t »« 


x 2 + r Q 2 - 2xr Q cos 0 O + t 


0 


/Xr 0 COS 2 expCik^+r 2 -2xr. cos 0. + £ 2 ) 

f ; ^—2 5 0 . . . . d £ 


+ 2^ r 0 Sln 6 0 


0 


(x 2 +r 0 2 -2xr 0 cos 0 Q + £ 2 ) 3 ^ 2 


where the first term alone exhibits a singularity at x = 0 . Hence the full solution 

A 

of (31), V(x) » V(x) + V(x) , remains finite in the limit x Of if 

• e ° 

. sm -r- ikr. . . .. 

— e - A(<+ k) e" 17r = 0 

/r o 

6 0 

sin ~y i(kr Q +7r/4) 


or 


v / rrr^ /<+k 


(37) 


This agrees with a determination by Jones utilized, in turn, for estimating the 
distant sound field amplitude along directions close to that of the wake and 
elsewhere. 

A different measure of the effect brought about by imposing a Kutta 


condition, namely the change in total energy radiated, can be given; in order 
to arrive at this, consider the integrals of an identity obtained from the first 
equation in (28) , 
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~2 I® V ' (<P V<J>) Im{6(x-x 0 )6(y-y 0 )(() } , 

which are taken, separately, over the adjacent domains y < 0 with a common 
boundary at the half-plane and the wake. After replacing the (two-dimensional) 
divergence integrals by curvilinear ones along great semi-circular arcs in y < 0 
and line integrals along y - 0+ , adding the results and employing the last 
two equations of (28), it turns out that 

00 * 

1 + 4 Im 4> s (x 0 ,y 0 ) + 4 Im / A e" 1KX V(x) dx (38) 

0 

where the (familiar) term involving the value of the secondary wave function 

<J> s (x,y) = <f>(x,y;P) - i (k/(x-x 0 ) 2 +(y-y 0 ) 2 ) 

at the fixed source point p ( x Q>yg) is supplemented by a non-local term, on 
account of the postulated wake or discontinuity line for the velocity potential. 

A 

There is no contribution to the latter term stemming from the function V(x) 
and, specifically, its Fourier transform (cf. (33)) 

V(-k) = f e _lKx V(x) dx = - ^ 

0 2 ie 

if K > k and >4 c-k has a wholly imaginary value, for then 

Im A V(-<) = 0 • 

This means that the Kutta condition affects terms in (38) which depend linearly 
on the complex magnitude A and appear in the combined form 
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A &= 4 Im |- | / H^ 1) (k/(x-x 0 ) 2 +y 0 2 )V(x) dx + A* / e" lKX V(x) dx 

2 f “ y 0 . ] 

Im < i / — V(?) d£ L + 4 Im A V(-k) 


77 


whose details are shaped by the representation (30) and the Hankel function integral 


H< x) (k/(*-* 0 ) 2 +y 0 2 ) - - 

—00 


L5(x-x 0 )+i>4 2 -; 2 |y 0 

77 


It suffices here, in view of the crude model, to merely estimate 
when kr Q and <r Q are very small; thus, relying on the explicit determination 

A 

(33) for V(£) the first term of (39) becomes 

r -iCx 0 +iv4 2 -s 2 y n ^ 


^lm^if^d=k / 


C + < + ie .AT? 


and, since the integral remains finite in the limit x q»7q 0 > this has an 


approximate value 


iySklmiA 7 - — 2 r - 1 

" \ .» ? + * + ^ X?F J 


= 2 Im A , kx 0 0 


On retaining just the initial and predominant term in the representation (36) 


of V(x) and noting the relation 

;Wi 77 e lk(! * r o ) V i ikr 0 » e i<k+?, v 


4 Im A*V(-k) 


- f to { 


.* , e o r dt 

A sin -y J 

0 /r (1+t) 


0 * 

= 2 sin — Im A , <r^ ■+■ 0 . 


it follows that 



- 22 - 


Accordingly, 

9 0 

A 9 * - 2 Im (A-A sin — ) 

0 Q 

sin • (1-sin -~) 

= , (<+k)r « 1 (40) 

*2 (<+k) r Q 

which result indicates that enforcement of a Kutta (or finite velocity) con- 
dition at the edge of a half-plane significantly elevates the power output from 
a fixed line source situated nearby, if 0 Q is not too small. 


§5 . Radiation by a Point Source Above a Plane Impedance Boundary 

Consider an infinite plane surface S which is locally reflecting and 
characterized by the uniform normal impedance pcC ; then, if z = 0 on the 



plane, the boundary condition 


+ ikr)4> ■ 0 , z = 0 (41) 


applies, with r) — — termed the specific 
admittance. It follows directly from (41) 
and the normal (z-) component of the flux 
vector (6) at the surface that 

N z Im(< = ~ Im(ikn|<f>| 2 ) 


and thus energy leaves, rather than enters 
the half-space z > 0 if 


Re ri > 0 . 


(42) 


Let the equation (4) describe a source at P(0,0,h > 0) and introduce 


the split 
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<t>(r) - <f> p (r) + <(> s (r) 
ik | r-P | 

■ T, + M?) » z > 0 (43) 

4 tt | r-P | 

where <f> s (r) is regular everywhere above the surface S . Multiplication of 
& 

(4) by <j> and subsequent integration over the half-space z > 0 leads to the 
relationship 

9 = 1 + ~ Im <|> S (P) = S+Zfi (44) 

in which 

$= 4tt Lim {r 2 / |<J>| 2 d^} (45) 

r-*» 

represents the energy outflow rate in all directions above the impedance sur- 
face and 

SJ = 4tt Re {r| / dx dy} (46) 

S 

represents the rate of energy loss or absorption at the complex impedance surface, 
both of these measures relative to the power output from an identical source in 
a completely homogeneous environment. 

If the component parts of the wave function (43) are given explicit 
and individual expression through the integral 


<f)(r) - <f>(R 


>z ) - _L f - V gR) I'g.-^-k 2 1 z-h | [ >4 2 ~k 2 + ikn ^-^ 2 -k 2 (z+h) 

n / o o I /n 1 a 


^Tk 2 




d? 


ikq 


arg 



0 , K > k 
-tt/ 2 , £ < k 


(47) 


that takes account of both the cylindrical symmetry round the line OP normal to 
the surface and the boundary condition (41) , it appears that 
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<t> q (P) = <p (0,h) / 


0 


2 -k 2 + ikn e -2v4 2 -k 2 h d? 


r-k - ikn 


and hence 


Im <t> (0,h) = § / ^- e - X P ( ~ 2 ^ 2 - k2 ^ - d£ + ^ / 
2 " k 2 "° 


M f £ sin(2^ 2 -g 2 h) d5 
0 (v4 2 -£ 2 +ka) 2 +(k6) 2 


kV + kV + g 2 - k 2 


0 yfi-E 2 (»4?-£ 2 +ka) 2 + (kB) 2 


( 2 v { 2 -£ 2 


£ h) d£ (48) 


where a,B denote the real and imaginary parts of the admittance, i.e., 


n = a + iB 


On referring to the pair of Fourier— Bessel integrals which stem from 
(47) , namely 

1 

<J> ( R»0) = - - / F(£)£J 0 (£R) d£ 

and 

00 -*^ 2 -k 2 h 

F(£) = / <J>(R»0)RJ (£R) dR = — - e (50) 

° v^ 2 -? - ik(a+±B) 

it is next found that 


/ |<f>(R, 0) | dx dy == 2 tt / | cf> (R, 0) | 2 R dR 

-°o 0 

°o oo 

=-i / RdR / ££ , J n (£R)J 0 (£'R)F(£)F*(£ t ) d£ d£' 

0 0 u u 

= /6(£-£’)£ F(5)F*(5‘) /£|F(£)| 2 d£ 


i / _ ^exp (-2^ -k h) d5 + 1 , gdL_ 

k (v^ 2 -k 2 +kB) 2 + (ka) 2 V 0 (y^Ip+ka) 2 
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All the integrals encountered in (48) and (51) , save for 


k 

I - / 




= ; 


ydy 


0 (^ 2 -S 2 +ka) 2 +(k8) 2 0 |n+/V| 2 


_ 1 (1+a) 2 + 8 2 a 

■ 2log "7T7""W 


(tan 


-1 1 + a 


w 


- tan 


-I a 


TeT ) 


are reducible to the exponential type, i.e., 

00 — £ 

E 1 (z) = f - - , |arg z| < tt ; (52) 

z ^ 

and the consequent determinations of the total power factor , radiation 
factor $ , and absorption factor gJ prove to be 

S’ - 1 + 3i ” k ^ kh + 2Re{t 1 e~ 21khn E 1 (-2ikh(l+Ti))} , (53) 


s - 1 + sa_2Si . 2ct / — Hit 


2kh 


|n + v / f-y 2 1 2 
+ 2 Re{ne i2lkhn (E 1 (-2ikh(l+n))-E 1 (-2ikhn))} , 


and 


(54) 


0=2a / ^ + 2Re{ne _2lkhn E (-2ikhn)} (55) 

0 |n + v'Cy 2 | 

with proper limiting values if a , 8 0 or a , 8 00 (the hard or soft surface 

conditions, respectively) . 

Ingard and Lamb state that "a rigorous solution is not so easily ob- 


tained (for absorbing surfaces) since the evaluation of the reflected field 
cannot be obtained directly by simple image procedure but would have to draw 
upon the analysis °f reflection of a spherical wave by a boundary.” Ingard* s 
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own study of the latter problem (1951) , which expresses the complete wave 
function as a sum of primary and "image" source contribution (see figure) 


*(r) 


ikr. 


4ttt- 


+ Q 


with a complex and variable amplitude factor 


Q(r 2 ,0) = 1 - 2kr 2 n / 


ikr„ 


4frr 2 


"kr C 
e f d£ 


( 56 ) 




ri cos 0 + i£) ^ + sin 2 0 (1— rj^) 


(57) 


affords, it may be noted, a ready determination of SP ; thus, from (44), 

2ikh ' 


& = 1 + 


X im|q(2h. 


0) 


4TT*2h \ 


, sin 2kh „ T I 2ikh . 

1 + — — 2 Im < ne / 


-2kh£ 


= 1 + 


2kh 


sin 2kh 
2kh 


1 + n + i? 


r? 


+ 2Re|ne“ 2lkhn E 1 (-2ikh(l+n))j 


in agreement with (53) . The wave function (56) is, on the other hand, less 
suited than the version (47) in terms of cylindrical coordinates (R,z) for 
determining the absorption factor . 

§ 6 * Power Output of a Point Source in a Uniform Flow 

Acoustical excitations and their power indices are strikingly affected 
by the existence of a background flow or a relative motion between the source 
and its environs. Assume that the convected wave equation 


V 2 $ " ~T + ^ ‘ V)2$ = -<5(r-P)f(t) 


(58) 
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characterizes the velocity potential for a monopole source at the fixed point 
*> 

P , with temporally variable amplitude rendered by the function f(t) , in a 
fluid which streams uniformly and steadily at velocity U . Multiplication of 
(58) by -p3$/3t , where p designates the constant density of the fluid in the 
absence of acoustical perturbation therein, and the subsequent recasting of 
individual terms leads to the relation 


where 


and 


fr+ V • N = p6(r-P)f(t) || 


W = \ p(V$) 2 + -fij (||) 2 - (U • V$) 2 
2 2c 2 9t 2c 2 


N = -pf v* + -^§ (f+U-V*) . 
c 


(59) 

(60) 
(61) 


The scalar and vector quantities W , N may be identified as measures of the 
acoustical energy density and flux, respectively, the latter being directly 
applicable to calculations of energy transfer across fixed surfaces in the fluid. 

When P lies at the origin, f(t) = e and the streaming velocity 

U : (U,0,0) is along the positive x-direction, the time-periodic (complex) 
source function deduced from (58) , 


$(r,0,t) 


ikR 
e 

4ttR^ 


-icot 

e 


k = u)/c 


(62) 


has both non-symmetric amplitude and phase factors with the specifications 

-Mx + R„ 


2 2 2 2 2 
R x - x + (1-M ) (y +z ) , R = 


1 - M 


(63) 


at Mach numbers M = U/c less than unity. 
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The determination of energy output from this source by integrating 
the flux vector N is facilitated with a suitable choice of the enveloping 
control surface; and an evident one comprises a cylindrical (curved) portion, 
aligned with the flow direction and having the source on its symmetry axis, 

/2 2 

together with plane end sections normal to the flow. If a - V y +z = constant 
defines such a cylindrical surface, the normal flux thereat is 


N 


a 


M 3$ 
P 9t 9a 


and thus, taking into account the representation 


M.lJ (1-M 2 )a _a_ 

9a - 4 tt 1 R 3 + lk 2 

R 1 R 1 


exp(-iu(t --)) 


the time-average measure 


N 


a 


1 / l x 2 pk^ a 

2 V c _ 3 


follows* Since a vanishingly small energy flux passes outwards across finite 
plane areas perpendicular to the cylinder axis at remote sites (where |x| -> 00 ) , 
the whole power delivered by the source has the magnitude 


& - 2tt / N a dx = TTp — (~) 2 
.00 a c 

= mi i 
8irc i - m 2 


2 J 

a dx 

r 2 /. 2. 2.3/2 

[x +(1-M )a ] 


(64) 


An identical result is secured on integration of the right-hand side, 
or source term, in the time-averaged version of the energy balance equation 
(59) , namely 
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<?-^Re / p<5(x)<$(y)6(z)e iaJt || dx dy dz 


| po) Im / 6(x) 6(y) 6(z)$(r,0,t)e^ lJt dx dy dz 


H im / exp 


ika \ <$(y)6(z) dy dz 


/-M 2 / X-M 2 


_ Lim exp^Cika/yl-M 2 ) _ pa) 2 1 

" 8,1 <«> ' ’ 8 ” C 1 - M 2 ' 


Consider next an identical source at distance h above a rigid wall 
occupying the entire plane z = 0 ; then the relevant wave function is composed 
of two like terms 


$ = $(r,P,t) + $(r,Q,t) , z > 0 


where P(0,0,h) and Q(0,0,-h) are mirror image points in the plane, and the 
radiated energy can be directly found by recourse to a source integral of the 
form (65), viz. 

00 

& = | po) Im / 6(x)6(y)6(z-h) ($(r,P,t) +$(r,Q, t))e itUt dx dy dz 
—00 

_ pu) 2 1 pw sin(2kh//l-M 2 ) 

" 8 ire T ~2 + 8 tt 7== * (67 > 

1 " M 2hA-M 2 

The relative power factor obtained from the ratio of (67) and (64) 


-1 + £I 
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which incorporates the effects of both a background flow and a hard reflecting 
surface, generalizes the earlier expression (3). 

It requires a rather more involved procedure to arrive at (67) by 
calculation in the far field; thus, introducing the spherical polar coordinates 
r,0,ip through the transformation formulas 


- r cos 0 , y = r sin 0 cos \p , z = r sin 0 sin ip , 


with the ranges 


0 < 6 < TT , 0 < \p < TT , r > 0 


appropriate to the half-space z > 0 , and employing the asymptotic character- 
izations 


R_ = v4 2 +(1-M 2 ) [y 2 +(z+h) 2 ] = /-M 2 Sin 2 0 r + h(l-M 2 ) — — - 9 . sirx . 1 

+ / 2 2 

»a-M sin 0 


r oo 


-Mx + R 


±- , —I— J(X-M 2 sin 2 0 - M cos 0)r + h(l-M 2 ) - . 9 sin i/j 


1 - M 2 1 - M 2 


an estimate for the complete wave function, 


A ~-~ 2 • 2 


•M sin 0 J 


* 1 cos(kh sin 0 sin ^/>q-M 2 sin 2 0) I., /1-M‘ sin“ 0-M cos 0 

w ' 2irr . . ex P <{ ik 9 r-iutj. , 


A -K 2 * 2 

fcJAp J. K. 

v 4-M 2 sin 2 0 ' 1 - 


r +«o (69) 


is obtained from (66) . On employing the latter along with the radial component 
of the energy flux 

N r - -p {(1-M 2 cos 2 6) |f + MislnjLcos_e |£ _ a cos e |£ } 


30 c 


3t 



-31- 


the average radiated power in the half-space z > 0 acquires the versions 


IT IT 

& = / dip f sin 0 d0 Lim (r N ) 
0 0 r-*» r 


= 7 a* / sin e d 6 ■! -fi (i) 2 ot>s 2 [ kh JiULaiaJt. 


0 


/-M 2 


sin 2 0 


[‘ 


Vl-M 2 sin 2 0 - M cos 0 1 - M 2 cos 2 0 


1 - M 


2 2 

1 - M sin 0 


+ k 


M sin 0 cos 0 1 


i M 2 , 2 n , J d 0 

1 - M sin 0 1 - M 


d (v^-M 2 sin 2 0 - M cos 0) 


+ kM 


cos 0 


2 2 
1 - M sin 0 


£<s> 2 / d * / 


sin 0 


n n ^2 . 2 m 3/2 

0 (l-M sin 0) 


J / 


2 / , , sin 0 sin ^ i 
cos I kh ^ |d0 


"sin^ 0 


TT 


JTp , 0 ) x 2 ' r ' sin 0 


4c 


/ w \ ^ r am 

' 2 tt J 2 


0 ( 1 -M sin 0 ) 


2 AN 3/2 *1 1 + J 0l 2kh 


sin 0 


sin 2 0 /J 


d 0 


(70) 


where denotes the zero order Bessel function. 

That part of (70) which does not involve the Bessel function agrees 
precisely with the first term in (67) ; to link the remaining parts a change of 
the integration variable proves convenient. If 0 is replaced by y in (70), 
such that 

sin 9 sin V _ 0 < 0 < n/2 

■4-M 2 sin 2 e vd 2 0 < Y < n/2 

and 


sin 0 = 


sin y 


cos 0 d 0 cos y dy . n / J2 . 

, 5 5 j /2 = ~~ > tan ® " v^L-M tan y , 


►4-m 2 cos 2 e (1 -^ sln 9 > 



then 
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/ 


sin 9 


0 (1-M 2 sin 2 6) 3/2 0 


tt/2 

= 2 / 


2kh 


sin 9 


/-M 2 sin 2 9 
sin 9 


d9 = 


0 (1-M 2 sin 2 9) 3/2 0 


J~ 2kh 


sin 9 




✓1-M 2 sin 2 9 I 


I d9 


tt/2 


2 / sin y J 2kh — — Y dy 

°\ vG? ' 


1 sin 


(2kh/X-M 2 ) 


kh 


since 


ir/2 

/ sin x J (a sin x) dx = 


sin a 
a 


and the full compatibility of (67), (70) is demonstrated. 
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